Preprint typeset in JHBP style - PAPER VERSION 



Flavour Structure of R-violating Neutralino Decays at 
the LHC 



o , 

CN , Nils-Erik Bomark, 

Department of Physics, University of Bergen, Postboks 7803, N-5020 Bergen, Norway 
E-mail: Nils-Erik. BomarkOift .uib.no 

CN ' Debajyoti Choudhury, 



Ph. 



Dept. of Physics and Astrophysics, University of Delhi, Delhi 110007, India 
E-mail: debajyoti . choudhurySgmail . com 



Smaragda Lola, 

CD ■ 

(~| ' TH Division, Physics Department, CERN, 1211 Geneva 23, Switzerland and 

Department of Physics, University of Patras, GR-26500 Patras, Greece 
^SJ ' E-mail: Magda.Lola@cern.ch 

> ■ 

Cn| ■ Per Osland 

(N ■ 

' Department of Physics, University of Bergen, Postboks 7803, N-5020 Bergen, Norway 

■ E-mail: Per.Dsland@ift.uib.no 

o 

Abstract: We study signatures of R-parity violation in the production of supersymmetric 
particles at the LHC, and the subsequent decay of the hghtest neutrahno being the end 
product of a supersymmetric cascade decay. In doing so, we pay particular attention to 
^ \ the possible flavour structure of the operators, and how one may discriminate between 

\ different possibilities. A neutralino LSP would couple to all quarks and leptons and a 

comparative study of its decays provides an optimal channel for the simultaneous study 
of all 45 R-violating operators. By studying the expected signals from all these operators, 
we demonstrate the ability to understand whether more than one coupling dominates, and 
to map the experimental signatures to operator hierarchies that can then be compared 
against theoretical models of flavour. Detailed comparisons with backgrounds, including 
those from MSSM cascade decays are made, using the PYTHIA event simulator. 
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1. Introduction 

Over the last decades, there has been increasing interest in R-violating super symmetry, as a 
plausible alternative to the Minimal Supersymmetric Standard Model (MSSM). Supersym- 
metrizing the Standard Model (SM) does allow for additional terms in the superpotential 
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over and above those corresponding to the Higgs potential or the Yukawa terms within the 
SM. The supersymmetric versions of the latter read 

IJ.H1H2 + meLiEjHi + mdQiDjHi + my,QiEjH2, (1.1) 

where Hi^2 are the Higgs superfields, L (Q) the left-handed lepton (quark) doublet super- 
fields, and E {D, U) the corresponding left-handed singlet fields. Gauge symmetry as well 
as supersymmetry do allow for additional bilinear terms of the form 

liiLiH2 (1.2) 

as well as trilinear ones, namely 

XijkLiLjEk + X'ijk^iQjDk + X'-jj^UiDjDk. (1.3) 

Each of the terms in eqs. (1.2) & (1.3) violate R-parity. In particular, the existence of 
any of l^'i, Xijk, X'-jj^ implies lepton number violation while a non-zero X'-^f. violates baryon 
number. SU(2) and SU(3) invariance imply that there are 48 R-violating couplings in all 
(3 /ij's, 27 A^^s and 9 each of Xijk and X'-jj.). 

The strictest bounds on lepton and baryon-number violating operators come from 
proton stability. The assumption of a conserved R-parity automatically rules out all of the 
terms of eqs. (1.2) & (1.3), rendering the proton stable [1], modulo higher dimensional terms 
endemic to the MSSM. However, alternative symmetries, namely baryon or lepton parities 
[2,3] can also exclude the simultaneous presence of dangerous LQD and UDD couplings [4]. 
Experimental constraints from the non-observation of modifications to Standard Model 
rates, or from possible exotic processes [5-9] also impose additional limits. Overall, the 
phenomenology to be expected out of such theories is very rich, since the LSP (Lightest 
Supersymmetric Particle) is no longer stable and the missing-energy signatures of the 
MSSM [10] are substituted by multi- lepton or multi-jet events [11-13]; single superparticle 
productions [14, 15] are also possible. 

In addition to the consequences for collider searches, R-violation implies that gravitinos 
(which may have been thermally produced after a period of inflation) are also unstable. 
However, gravitino dark matter in the framework of R-violating supersymmetry is plausible 
[16-18], provided that the gravitino decays slowly enough for its lifetime to be larger than 
the age of the universe [19,20]. This is an exciting possibility that allows for supersymmetric 
dark matter, even if the R-violating couplings axe sufficiently large to lead to observable 
signatures at colliders [17,18]. Moreover, it was found that the branching ratios for gravitino 
decays are very sensitive to the flavour structure of R-violating operators [17, 18]. 

In general, the flavour structure of R-violating couplings is of particular relevance in 
defining the nature of the signals to be expected and any information on it would be crucial 
for understanding the flavour structure of the fundamental theory. Indeed, one may try to 
relate hierarchies amongst R-violating couplings to those in fermion masses [21,22], using 
models with family symmetries. For example, a large class of such models allow only the 
third generation fermions to be massive, while the remaining masses are generated by the 
spontaneous breaking of the family symmetry. In such a scheme, if R parity is violated. 
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couplings with different family charges are hkely to appear with different powers of the 
family symmetry-breaking parameter, and thus with different magnitudes. 

Testing the above models against observations is very hard, since by allowing the most 
general family structure of couplings one ends up with a large number of possibilities. Single 
superparticle productions can lead to spectacular signatures, but they are very dependent 
on the flavour structure and one must ensure that the initial scattering states couple to 
the operators of interest in a specific model. If this is not the case, the respective single 
superparticle production mechanisms will be invisible. 

Prom this point of view, pair production of superpartners and their subsequent cascade 
decays via an unstable neutralino LSP have a great advantage, since the latter (by coupling 
to all quarks and leptons) could decay via any of the 45 trilinear operators, thereby allow- 
ing a comparative study. Through a detailed, correlated study of these decay chains, one 
may also investigate whether more than one R-violating couplings are of substantial size, 
"map" their magnitudes and hierarchies, and compare against theoretical models. This is 
the aim of the present paper, which is structured as follows: in Section 2, we summarise 
the generic framework of the analysis. In Section 3, we study energy and invariant mass 
distributions, for LLE operators, as well as leptonic R-violating branching fractions. In 
Sections 4 and 5 we discuss energy and invariant mass distributions for LQD and UDD op- 
erators, respectively. In section 6 we focus on the correlations between different couplings 
and possible mappings between the kinematic distributions and operator flavour hierar- 
chies. In section 7 we present our conclusions. Finally, in Appendices A and B we present 
spectra from three-body decays and theoretical invariant mass distributions, respectively, 
integrating over the unobserved neutrino energy. 

2. Framework for the Analysis 

We study possible effects of R-parity violation (RPV) in the end stages of cascade decays 
of supersymmetric particles. In other words, we assume that the RPV-couplings are small 
enough not to materially change either the production processes (which continue to be 
gauge-interaction driven), nor in the cascade decays down to the neutralino LSP. The LSP, 
though, decays, within the detector, on account of a non-zero RPV-coupling. A neutralino 
LSP is motivated even in the case of R-violating supcrsymmctry, since, due to the lack 
of electromagnetic or colour interactions, its mass after renormalisation group runs tends 
to be smaller than those of other sparticles for a very broad region of the supers3mimetric 
parameter space. A sneutrino or a stau LSP (due to the large Yukawa coupling, staus can 
be anomalously light) is still plausible, especially if mo is sufficiently small compared to 
mi/2- Those scenarios are beyond the scope of this work. 

Depending on the scenario, a variety of supersymmetric particles could be produced 
in proton-proton collisions, mostly squarks and gluinos, 

qq qq, qg qg, gg gg, gg qq, qq g|, (2.1) 

but also charginos and neutralinos. The squarks and gluinos will typically decay to quarks 
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Figure 1: Feynman diagram describing one important supersymmetric decay chain followed by an 
R-parity violating neutralino decay. 

and leptons, one important channel being 



The importance of this is due to the presence of leptons, which allows easier extraction 
from the background. However, when confronted with leptonic or semi-leptonic R-violating 
neutralino decay, such leptons will constitute a background for our signals. 

We shall, in turn, allow the LSP to decay via LLE, LQD and UDD couplings, as 
depicted in Fig. 1. 

Our analysis has been performed for some of the SPS points [23]. These points have 
been identified so as to satisfy the WMAP constraint on dark matter [24] . While we do not 
require stable dark matter, these benchmark points are, nevertheless, convenient reference 
points in the multidimensional parameter space and allow for direct comparisons between 
the MSSM results in the literature and the expectations in the presence of R-violating 
operators. Several of the results we obtain are similar for all these points and therefore, 
for many of the distributions to be discussed, we focus on the SPSla case. However, for 
those predictions that are more sensitive to the SUSY parameter space, we present more 
global results, elaborating on the differences to be expected in each case. 

In the above cascade decay, the actual strength of the R-violating coupling does not 
lead to qualitative differences in the predictions, as long as: 

(i) it is strong enough for the neutralino decays to be inside the detector (> 10~^ for 
100 GeV sparticles and scaled accordingly for higher masses [25], in the case that we 
do not have any additional effects from phase space suppressions); 

(ii) it is not sufficiently large for sparticles to decay directly via the R-violating operator 
instead of reaching the end of the chain and then decaying via neutralinos. 

For couplings closer to 10~^ there might be displaced vertices that can improve the extrac- 
tion of the signal [26,27]. This, however, has not been studied in this work. 




qqt 




(2.2) 
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The numerical value of the upper limit of couplings does depend on the SUSY param- 
eter space and the flavour of the R-violating operator under consideration. In Section 2.1, 
we show that the cascade decay actually dominates for a large part of parameter-space 
as long as the R-violating coupling is kept within experimental bounds. Especially when 
more than one coupling is large, the strong bounds imposed on products of couplings from 
flavour violating processes [28] will assure cascade dominance. 

To analyse the above scenario, we simulate proton-proton collisions at a center of mass 
energy of 14 TeV using the PYTHIA 8.1 [29] Monte Carlo event generator. In order to 
correctly describe the kinematical distributions, we have modified PYTHIA to use the 
energy distributions of appendix A for the neutralino decay. For the SUSY input we use 
the SUSY Les Houches Accord (SLHA) [30] spectra from SoftSUSY [31], together with 
SLHA decay tables calculated in PYTHIA 6.4 [32]. 

The jet analysis is done using Fast Jet 2.4 [33] with a fcf-algorithm [34] where the param- 
eter i?jet) denoting the upper limit on AR = {A-qY + (^^)^ for a merge to be allowed, 
is set to 0.4 unless stated otherwise. Included in the analysis are all detectable particles, 
i.e. all particles except neutrinos (and stable neutralinos and gravitinos if applicable). 

For analyses including taus and b-quarks, the definition of tau (b-) jets is that at least 
60% of the pt of the jet should come from decay products of the tau (b-quark). This 
fraction is known in the Monte Carlo but for a more realistic study one would need to 
resort to proper tau and b tagging. For the rest of the paper we shall denote the tau-jets, 
Tjet and the b-jets, 6jet, except in figures where we use r and h in order to save space. 

A lepton is considered isolated if there is no jet with px > 10 GeV within a cone of 
AR < 0.4 and the total additional energy within AR < 0.2 is less than 10 GeV. Note 
that leptons are included in the jet clustering, but jets that, after the clustering, contain 
essentially nothing but a lepton, are not included in the further calculations involving jets. 

We also impose acceptance cuts, requiring pT > 5 GeV for leptons and pT > 10 GeV 
for Tjet's and bjet's- Where light quark jets are used, they are required to have pr > 20 
GeV. The reason for the harder cut on light jets as compared to rjet's and 6jet's, is to 
remove some background jets from beam remnants, bremsstrahlung etc. 

All data presented is based on simulations of one million events for each operator/SPS 
point combination. The backgrounds, tt, W, Z and QCD (with total px > 300 GeV) are 
studied simulating 10 million events for each in order to account for their larger cross- 
sections. 

2.1 Cascade dominance 

In order for the cascade to proceed to the end, giving an (unstable) LSP, the competing R- 
violating decay of the squark cannot be too strong. The decay of the squark to a quark and 

a neutralino (or chargino) proceeds via a combination of gauge and Yukawa couplings, in 
competition with an R-violating decay to a quark and a lepton with a strength proportional 
to some A' squared^. The relative rates of these channels will depend on the squark mass, 

^Analogous statements hold for decays through A", and in case of slepton/sneutrino decays for A- 
couplings. 
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Figure 2: Dominance of MSSM decay of squark vs its R-parity- violating decay, shown for two 
squark masses, and two values each of tan /? and A'. 



as well as on the neutralino and chargino spectrum, and the coupling A'. In order to 
illustrate this competition, we determine neutralino and chargino spectra in terms of M2, 
H and tan/3 (holding Mi = (5/3) tan^ 0i4/M2) and plot, in Fig. 2, the ratio 



Er(g^gx.") + Er(g^ q'x 

T{q^ all) 



(2.3) 



Where this ratio is high (represented in red and yellow) , the decay of the squark will proceed 
via a cascade to the lightest neutralino, which then will decay via the R-parity-violating 
interaction. Note that we do not need any assumption here about where in the decay chain 
the studied squark is produced and hence a high ^ for all sparticles present in the coupling, 
will ensure that the cascade continues down to the lightest neutralino. 
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The white regions in these panels are excluded by the lightest chargino mass being in 
conflict with the LEP bound, M-± > 94 GcV [35]. 

For some representative combinations of the squark mass M, tan/? and A', we see in 
Fig. 2 that even for rather substantial R-violating couplings, the cascade decay dominates 
for large regions of parameter space. One should also remember that typical decay chains 
may involve particles other than those present in the dominant R-violating coupling and 
thus have no open channels for direct R-violating decay. Since these cascades will always 
proceed to the lightest neutralino, the probability of interrupting the cascade with a direct 
decay, is further reduced compared to the above discussion. 

3. LLE couplings 

We shall first consider the couplings LiLjE^, which, typically, lead to decays of the LSP 
to two charged leptons and a neutrino, for example 

Xi — t e+iJ.~iye, e~/i+Pe, e'^e'Uf^, e^e~v^. (3.1) 

While isolated Icptons in the final state would also arise in R-conserving scenarios from 
decays higher up in the cascade chain, we will see that differences in kinematic distributions 
allow us to statistically disentangle the composition of the final state. 

We assume that only one LLE operator is non-vanishing and set the corresponding 
coupling to be 

Xijk = 10-'. (3.2) 

It should be noted that, even if two (or more) couplings of such a size are introduced, the 
constraints from FCNC [28] are satisfied almost trivially. There do exist a few combinations 
for which even such sizes can be of concern, but these we take care of specifically. Note also 
that many such constraints actually apply in combinations like ^ijk^iin/^i^ ^^"^ given 
that the couplings A can assume either sign (and, in principle, any phase), larger couplings 
could, in principle, be accommodated if one were to admit cancellations. Collider signals, 
on the other hand, would be insensitive to such cancellations. We, however, desist from 
considering such possibilities. Identical observations would hold for A' couplings as well, 
while, for A" couplings, the bounds are, in general, even weaker. However, wc will return 
to the issue of a comparative study of the magnitude of the operators, if several of them 
are large — a main aim of the paper — in a subsequent section. 

For the analysis of the LLE operators, the following selection cuts have been used [36]: 

• at least three isolated leptons with > 20 GeV; 

• the leading lepton has pt > 70 GeV; 

• pf'' > 100 GeV. 

These cuts are sufficient to reduce tt events by a factor ~ 3- 10~^ and W and events 
(including all production of weak bosons) by a factor « 6 • 10~^. The RPV events, on the 
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other hand, are reduced to approximately 1 — 40% of the total, with the suppression factor 
depending on the point in the SUSY parameter space as well as the nature of the RPV 
interaction. For operators with a lot of tau flavour, the reduction is significantly stronger 
than for operators involving only light flavours. 

This means that in the case of a SUSY parameter point with low cross section and 
low lepton production in combination with a tau-rich LLE operator, wc might get some 
W and events in the sample, but detailed simulations show that this contamination in 
the studied distributions is not a significant problem. 

Overall, in the presence of R-violation via LLE operators, we can see spectacular 
signatures (final-state flavours, like-sign leptons, etc.) [13]; these become more pronounced 
for the L2L3E1 and L1L3E2 couplings, which mix all flavours. 

3.1 Final-state flavour correlations 

Since we cannot detect the neutrinos or their flavours, the characteristics of various decay 
channels have to be described in terms of the charged leptons in the decay. From each 
neutralino we get one pair of charged leptons and their identities are correlated with the 
flavours in the operator. 

It is important to remember that every oper- 
ator gives rise to two distinct channels, a conse- 
quence of the fact that one of the charged leptons 
is associated with the E superfield and the other 
with the L superfield. Since SU (2) invariance re- 
quires the L-type fields to have different flavours, 
this leads to exactly two different channels, mod- 
ulo charge conjugation. The operator L1L2E1, 
for example, will give both e+e" and e/x pairs 
with equal frequency. 

A further complication arises when tau flelds 
are involved, then light-lepton pairs may arise Table 1: Possible final states for LiE cou- 
from leptonic tau decays. And since our final plings. 
states require isolated light leptons, the leptonic 

tau decays will be favoured over hadronic ones. Similarly, decay channels wherein the super- 
particles decay directly into light leptons would have a higher detection efficiency compared 
to those involving intermediate tau states. Thus, the composition of final states measured 
is not the same as given directly by the branching fractions but depends nontrivially on 
the supersymmetric spectrum as well as the flavour structure of the coupling. 
The open channels for each coupling are given in Table 1. 

3.2 pt distributions 

We first look at the pT spectra of the electrons (henceforth, "electrons" will refer collectively 
to both electrons and positrons) and muons as well as hadronic taus. We shall assume that 
hadronic taus (rjet) can be identified (with some efficiency). Note that in the rjets the 
neutrino is not included. 
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Figure 3: Electron, muon and rjot spectra for L1L2E1 (left panel) and L2L3E2 (right panel) at 
SPSla. Notice the dips in the spectra at 20 GeV; this comes from the fact that a lepton with 
PT < 20 GeV can only pass the cuts if there are at least 3 other high pr leptons in the event. One 
can also see features at 70 GeV due to the cut. 



The electron, muon and rjet spectra for L1L2E1 and L2L3E2 at SPSla are shown in 
Fig. 3. The two panels look the same apart from different relative normalizations of indi- 
vidual curves. These different normalizations can give insight into the flavour composition 
of the operator. However, since we do not know how much the R-conserving decay chain 
contributes and also since the overall number of leptons is more or less dictated by the cuts, 
the most interesting information here might be the difference in normalization between the 
electron and the muon spectra. Since the R-conserving chain is not expected to give any 
difference here, it might give some hint on the dominant flavour in the coupling. 

3.3 Invariant mass distributions 

One of the best tools for studying leptonic decays of the neutralino are invariant mass dis- 
tributions. The charged leptons from the decay will give rise to invariant mass distributions 
as calculated in Appendix B. Since the flavours of the involved leptons are determined by 
the flavours in the coupling, the channels in which these distributions are to be expected, 
can be deduced from Table 1. 

Let us illustrate how this works: L2L3E1 gives rise to e^^^ and e^r^ final states 
(see Table 1). Therefore, we expect a prompt contribution of the type (B.2) in the e^i 
invariant mass distribution as well as contributions of the type (B.ll) (from leptonic tau 
decays) in the e// as well as ee distributions. In addition, the e^r^ final state should give 
a contribution of type (B.12) to the erjet distribution. 

Since our goal is to determine the flavours of the coupling, we want to study all dilepton 
distributions available, as well as distributions including rjet's. We then want to use the 
presence of the various signals together with Table 1 to determine the operator flavours. 

Since we cannot tell which leptons come from the same vertex, we will have a problem 
with unrelated lepton pairs, i.e., pairs of leptons from different parts of the event. This 
combinatorial background will, in some cases, make it difficult to see the signal. However, 
since these pairs are not related, they will equally often have the same charge as the 
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opposite, and since we do not expect any signals in the same-sign distributions, subtracting 
those distributions turns out to be a very effective way of removing these backgrounds. 

The effect of the same-sign subtraction can be seen in Fig. 4 where we show the 
di-lepton distributions for L2L3E1 at SPSla both without and with same-sign (SS) sub- 
traction. Prom this we can see that the same-sign subtraction is especially important in 
order to extract the smaller signals where leptonically decaying taus are involved. Also for 
the distributions involving rjet's it is crucial to apply same-sign subtraction in order to see 
the signals. 

3.4 Backgrounds and uncertainties 

The interpretation of any detected excess as a signal is possible only if we understand well 
the (backgrounds from) R-conserving decay chains and also the effect of the cuts on the 
shapes of the distributions. The latter issue is expected to affect the accuracy of the fit of the 
theoretical distributions to the data, which makes, e.g., a reconstruction of the neutralino 
mass more challenging. The first issue, on the other hand, depends on the SUSY model and 
is especially important for the distributions including rjet's where the contributions from the 
R-conserving chains look very similar to the contributions from the neutralino decay. This 
can be seen in Fig. 4 where both the /urjct and the TjctTjct distributions show unexpected 
features, closely resembling the expected signals. On the other hand, the corresponding 
triangular distributions in the ee and distributions are more distinct and can therefore 
be removed from the analysis (practically invisible in Fig. 4) . 

There are, of course, backgrounds from standard-model processes too, but these are 
well handled by the event selection criteria. Other experimental issues include photon 
conversion as well as photon and lepton misidentification. The necessity of same-sign 
subtraction also implies a sensitivity to the ability to correctly measure the sign of the 
leptons, especially the corresponding charge of rjet's. 

For some choices of SUSY parameters, a real may be produced in the decay chain 
and this would show up in the same-flavour (opposite-sign) dilepton invariant mass distri- 
butions. However, since this only gives a peak at the mass, this effect can be easily 
accounted for. 

3.5 Comparison of Monte Carlo results with theory 

We now compare Monte Carlo simulations with the theoretical distributions derived in 
Appendix B. In what follows, the neutralino mass and the normalizations of the spectra 
are fitted to the Monte Carlo distributions, thereby allowing for a determination of the 
neutralino mass. 

As can be seen in e.g. Fig. 20 in appendix B, most of the expected distributions 
look very similar, the main differences being different kinematical cutoffs. This means we 
have to be careful with the fitting procedure to avoid too large degeneracies in the fitted 
parameters; e.g. a distribution of type (B.2) may look very similar to a distribution of type 
(B.ll) with a larger neutralino mass. 

To avoid problems with this, each distribution has only been fitted using the distribu- 
tions that are relevant due to the assumed operator flavours. For the relevant distributions 



-10- 




20 40 60 80 100 120 140 160 180 200 
M[GeVl 



0.02 
0.015 

0.01 
0.005 




SPSIa 






:L2i-3Ei 










- - ez 


r ^ 






/ \ 

: / \ 




-.-HI 


/ \ 






; ; \ 






'- 1 ^ \ 






- / .-.-'i-. \ 

" >'* Vn. 

- i" 












■ f 




1 , . , 1 


. . 1 . . . 1 . . . 1 . . . 



20 40 60 80 100 120 140 160 180 200 
M[GeV] 



Figure 4: Di-lepton invariant masses for L2L3E1 at SPSla, where the neutrahno mass is 97.0 GeV. 
Thin curves: without same-sign subtraction. Heavy curves: M{££') distributions after same-sign 
subtraction. Here, "r" refers to a hadronic tau. 



we also include a triangular distribution where the cutoff is assumed known (only for SUSY 
models where such a contribution from the chain exists), and a peak at the mass. The 
parameters to be fitted are then the neutralino mass as well as the normalizations of the 
various components. 

In a realistic scenario where the operator is not known but to be determined, the 
relevant distributions can be determined by a comparison with Monte Carlo data. In cases 
where this is difficult one can exploit the expected correlations between the distributions; 
if one distribution has a clear cutoff, the others cannot have neutralino masses significantly 
lower than that cutoff. Note that at least one of the distributions will always have a cutoff 
close to the neutralino mass. 

Assuming L2L3E1 to be the dominant R-violating term, various dilepton mass dis- 
tributions corresponding to the benchmark points SPSla, SPSlb and SPS6 are shown in 
Fig 5. The left panels show the pure dilepton invariant mass distributions while the right 
panels show the distributions involving rjct's. As one can see, the fits are in all cases quite 
good, but there are some deviations at the lower end, especially in the rjet distributions. 
This is due to the cuts imposed and the consequent deterioration of the signal profile. Es- 
pecially in the case of SPS6 we find that the inclusion of the triangular contribution from 
the chain into the fit is crucial for a correct fit to the smaller signals. 

The situation changes somewhat when the operators give large numbers of taus in the 
decay. However, as Fig. 6 shows, the fit is still rather good. Of special interest might be 
the case of L1L3E3 at SPS6 where there is no tau from the cascade decay chain and we 
can therefore see all the distributions expected, including the softer //Tjet signal where the 
fj, comes from a decaying tau. When comparing the various couplings at SPSlb it is clear 
that one cannot tell whether a signal in the rjetTjet distribution is from the neutralino decay 
or the decay chain, the expected distributions are much too similar. 

Table 2 shows the neutralino masses estimated from the fits shown in Figs. 5 and 6. 
Included are also statistical errors on the fitted neutralino masses. In addition we expect 
systematic uncertainties from the theory, especially concerning the inclusion of the event 
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Figure 5: Comparison between fitted theoretical distributions (thin curves in violet) and Monte 
Carlo results (heavy curves). The results for L2L3E1 are shown for three different SPS points. 



selection. Therefore, the uncertainties given in Table 2 should not be taken too seriously 
but rather as lower bounds on the uncertainty. 

It seems that the measurement of the neutralino mass from the erjet mass always gives 
an overestimate. The same is true for the di-lepton mass distributions when one of the 
leptons comes from a decaying tau, although, in the latter case, the overestimate is less 
severe. In the relevant cases, the distribution when both leptons come from decaying taus 
are included but the impact on the resulting fit is not very large; we see from the //// 
distribution for the operator, where both muons have to come from decaying taus, 
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Figure 6: The same as Fig. 5 for some tau-rich operators at SPSIb and SPS6. 



that the signal is too smah to get a good fit and the resulting neutralino mass estimate is 
completely useless. This is reflected in the very large statistical uncertainty. In general the 
possible improvement by including this last distribution is smaller than the uncertainties 
from e.g. the effects of the event selection. 

When interpreting these neutralino mass estimates, one has to keep in mind the large 
uncertainty associated with the effects of the event selection on the distribution. In the 
theoretical curves, these effects are included through the lower limit on the mass ratio 
entering in Eq. (B.ll). The value of this limit has great impact on the resulting distri- 
bution and therefore on the estimated neutralino mass. It is therefore clear that a better 
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model 


xY 


ee 




e/i 


er 


UT 


TT 


L2L3E1, SPSla 


97.0 


116 ± 0.3 




96.8 ± 0.1 


106 ± 0.5 






L2L3E1, SPSlb 


162.2 


167 ± 0.5 




171 ± 0.1 


171 ± 0.4 






L2L3E1, SPS6 


189.8 


195 ± 1 




187 ± 0.1 


200 ± 0.5 






L1L3E3, SPS6 


189.8 


205 ± 0.2 


295 ± 6 


195 ± 1 


198 ± 0.8 


227 ± 2 


224 ± 3 


L1L3E3, SPSlb 


162.2 


172 ± 2 


272 ± 13 


176 ± 0.4 


172 ± 0.7 


204 ± 3 


187 ± 3 


L1L2E3, SPSlb 


162.2 


173 ± 1 


172 ± 1 


172 ± 0.5 


172 ± 0.4 


172 ± 0.4 





Table 2: Fitted neutralino mass (all values are in GeV) from the various distributions. 

understanding of the effect of the cuts is very important in order to get the best possible 
neutralino mass estimate. 

In order to improve these fits one could make use of various expected correlations; e.g. 
as mentioned before, the neutralino mass of course has to be the same for all distributions 
and one could use the most suitable distribution, presumably the one with a cutoff at the 
neutralino mass, to constrain the fits to the smaller distributions with cutoffs significantly 
below the neutralino mass. Ideally one would want to use the fitted parameters to determine 
the decay channels for the neutralino, but due to the aforementioned degeneracy problems 
it is hard to achieve the required accuracy in the parameters, this could be improved if one 
uses the fact that e.g. an er decay channel for the neutralino gives correlated signals in both 
the erjet, ee and e/j, distributions, where the relative sizes of the signals are basically given 
by the branching fractions of the tau (together with tau tagging efficiency and detector 
acceptances). 

3.6 Some counting observables 

We now consider observables that can bring additional insight into the operator structure. 
One obvious candidate is the determination of the branching fractions of the neutralino 
into various channels. This can be achieved by studying the leptonic structure of the final 
states. 

To get a clean sample to work with, we look at the subset of all events where we have 
exactly four leptons of which two arc positively charged and two are negatively charged. 
This restriction increases the probability that the leptons have all come from decaying 
neutralinos. 

Let P2e-2e+ denote the fraction of these events that have exactly two electrons and 
two positrons: 

iV[(2e-)&(2e+)] 



2e-2e+ 



(3.3) 



(and similarly for -P2//-2/i+)) with 

N^i = 7V[(4e)] + iV[(3e)&(l//)] + iV[(2e)&(2//)] + 7V[(le)&(3/x)] + A^[(4/x)], (3.4) 

subject to the constraint that the charges should sum to zero. Further, let P2e2^ denote 
the fraction that has two same-sign electrons and two same-sign muons, i.e., 2e~2//"'" or 
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2e+2/Lt-: 

iV[(2e-)fc(2/x+)]+jV[(2e+)fc(2/x-)] 

In the event that all four leptons come from decaying neutralinos, Pe-e+ is the proba- 
bihty that a neutrahno decays to an electron-positron pair. Since under these assumptions, 
the only possibility to achieve 2e~2e''" in one event is to have two neutralinos that both 
decay to electron-positron pairs, we get 

P2e-2e+=Pe-e+- M 

Note that -Pe-e+ is a probability^ referring to one neutralino, whereas -F*2e-2e+ refers to the 
whole event. 

Similarly, ^2^-2^+ = P^-^+ and P2e2/. = ^Pef, where Pg^ = Pe-^+ + Pe+^- and 
Pjjt- fji+ iPe- n+ and Pe+ix- are defined in accordance with Pe-e+- The factor ^ in the second 
of these equations comes from the requirement that the same-flavour pairs also have the 
same charge. From this we can calculate the values Pe- e+ ^ Pfi- fi+ and Pg^. 

A problem with the above calculation is that it does not tell us anything about channels 
including taus. Since we require exactly four leptons, events with taus are only included 
when the taus decay leptonically. This also means that decay channels with light leptons 
are favoured by the cuts over those with taus and therefore care has to be taken when 
interpreting the numbers. However, keeping this in mind, the probabilities Pe- e+ ^ P/x- 
and Pe^ do tell us which decay channels are important. 

In order to get some information about decay channels involving taus, we need to 
also study events with three rather than four leptons. Let us study events with one pair 
of same-sign, same-flavour leptons {i^£^) and one extra lepton of opposite charge and 
different flavour {i'^). Such events must all contain one different-flavour, opposite-sign 
lepton pair (i^i'"^) and one extra lepton (i^), the extra lepton stemming from a neutralino 
that decayed cither to a r plus a lepton where the r decayed hadronically, or to a lepton 
pair where one of the leptons was lost to the detector. 

Since counting observables depend not only on production rate and detector efficiencies, 
but also on a priori unknown factors such as selection efHciencies, it is useful to study, 
instead, ratios such as 

P2,j.e N2^-e+ + N2fj,+ e- -^2/i-e+ ' 

where the numbers A'^ refer to the number of events with the specified number of leptons. 
The last equality above is based on the assumption that both charges are produced with 
the same probability and detected with the same efficiency. 

Let Tis now consider a particular case wherein we start with four leptons (e, //) stemming 
directly from the decays of two neutralinos with one of the four (randomly picked) getting 



^This is the probability that a randomly picked neutralino in an event with exactly four leptons, is 

decaying through a certain channel. In other words; it is a branching fraction, but depends strongly on the 
cuts used to calculate it and it should not be confused with a "true" branching fraction. 
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lost (absorbed in a jet, too little px, too large a rapidity etc.). The flavour distribution of 
the four-lepton state is then given by the probabilities P^-^+,Pe/i and -Pe-e+- Assuming 
that the efficiencies for muons are the same as for electrons^, we have 



where the factor 2 is a combinatorial factor arising from the fact that the individual leptons 
could have arisen from the two different neutralinos. 

Since the prediction of Eq. (3.8) does not take into account the possibility of the fourth 
missing particle being a hadronically decaying tau, it will, occasionally, give an inaccurate 
answer. This apparent lacuna could, however, be used to our advantage: a comparison 
of the prediction of (3.8) with a direct measurement of the ratio (3.7) would give us a 
quantitative measure of the tau decay channel. If the decay channel x5~^6 + ''" + ^isa 
significant one, the measured fraction will be much larger than the prediction (3.8), while 
if the channel Xi ^ + 7" + is present, Eq. (3.8) will give an overestimate. 

Independent measures of the branching fractions to various final states are given by 
the integrals over invariant mass, namely 



where 1,1' & {e, //} and the distributions are evaluated after same-sign subtraction. Whereas 
the P's above refer to the four-lepton events where the charges sum to zero, these integrals 
refer to all events that pass the cuts. 

Both the integrals over the invariant mass distributions and the probabilities Pe-e+ > ^11-11+ 
and Pe^ must be interpreted with the used cuts in mind; in both cases leptons are favoured 
over taus. The resulting branching fractions include Icptonically decaying taus in addi- 
tion to prompt leptons but the former contribution is strongly dependent on the cuts and 
therefore hard to estimate analytically. Due to the requirement of exactly four leptons, 
the probabilities -Pe-e+ -^e/* essentially include only events where there are no 
taus in the decay or all taus decay leptonically. The invariant mass distributions, on the 
other hand, arc expected to contain a larger amount of Icpton pairs from sources other 
than neutralino decay. This can, however, be handled by the same-sign subtraction as 
well as by comparing the distributions with the theoretical expectations. The probabilities 
Pe-e+j-f)* ^^"^ ^re more sensitive to contamination by single leptons from other 
sources and are therefore more dependent on a clean event sample. 

All the observables discussed in this section for all LLE couplings at SPSla are shown 
in Table 3. When comparing the probabilities Pe- e+ j n+ Pg^ and the integrals over 
invariant mass distributions, we see that they show similar structures, implying that they 
both are good indications for the relevant decay channels. 

In table 4, we compare these observables for the same operator {L2L^Ei) but at three 
different parameter points: SPSla, SPSlb and SPS6. As the values are similar at all three 
points, we may conclude that they are largely independent of the underlying SUSY model. 

^While this is not exactly true, the differences are marginal, and, at the level of sophistication of this 
study, can be ignored. 



-P2eM _ -Pe/^(2Pe-e+ + Pefx) _ 2Pe-e+ + Pe 



(3.8) 



P2^e Pe^.{2P^-^+ + Pe,.) 2P^-^+ + P, 




(3.9) 
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model 


Pe-e+ 


Pejj, 




(3.7) 


(3.8) 


J A/ee 








SPSla 


0.48 


0.51 


0.02 


2.5 


2.69 


0.48 


0.51 


0.01 


L1L2E2, 


SPSla 


0.02 


0.49 


0.50 


0.37 


0.35 


0.01 


0.49 


0.50 


T T IT' 


SPbla 


0.26 


0.48 


0.27 


0.93 


0.98 


0.27 


0.47 


0.26 




SPSla 


O.ol 


U.lo 


U.Ul 


1 n 
19. 


1 n 1 
iU.l 


O.OD 


n 1 
U.iz 


n no 
O.Oo 


LiL^E2, 


SPSla 


0.02 


0.81 


0.17 


0.26 


0.75 


0.02 


0.87 


0.11 


LiL^E^, 


SPSla 


0.45 


0.45 


0.09 


3.77 


2.15 


0.43 


0.44 


0.13 


L2L^Ei, 


SPSla 


0.17 


0.81 


0.03 


3.26 


1.33 


0.11 


0.87 


0.02 


L2-L3-E2, 


SPSla 


0.01 


0.16 


0.82 


0.05 


0.10 


0.02 


0.11 


0.86 


L2L3E3, 


SPSla 


0.09 


0.45 


0.45 


0.27 


0.47 


0.12 


0.44 


0.44 



Table 3: The three first columns give calculated values of Pg-g+, P^-^+ and Pe^ for all the different 
operators at SPSla. We also give the Monte Carlo value of the fraction (3.7) as well as the prediction 
from Eq. (3.8). Finally, we tabulate the integrals of the invariant mass distributions, normalized so 
that the sum is one. This should provide an independent estimate of the branching fractions of the 
first three columns. 



model 




p 




(3.7) 


(3.8) 




IK, 


IK, 


L2L3E1, SPSla 


0.17 


0.81 


0.03 


3.26 


1.33 


0.11 


0.87 


0.02 


L2L3E1, SPSlb 


0.18 


0.80 


0.02 


3.72 


1.38 


0.12 


0.87 


0.01 


L2L3E1, SPS6 


0.17 


0.81 


0.03 


3.54 


1.33 


0.12 


0.85 


0.03 


L1L2E3, SPSla 


0.26 


0.48 


0.27 


0.93 


0.98 


0.27 


0.47 


0.26 


L1L2E3, SPSlb 


0.25 


0.50 


0.25 


1.00 


0.99 


0.25 


0.50 


0.26 


L1L3E3, SPSla 


0.45 


0.45 


0.09 


3.77 


2.15 


0.43 


0.44 


0.13 


L1L3E3, SPSlb 


0.43 


0.49 


0.08 


3.37 


2.08 


0.41 


0.50 


0.08 



Table 4: Similar to table 3 for L2L3E1 at SPSla, SPSlb and SPS6, as weU as L1L2-E3 and L1L3E3 
at SPSla and SPSlb. 

However, also shown in table 4 are L1L2E3 and L1L3E3 at SPSla and SPSlb and here 
some noticeable differences arise. First of all, for L1L2E3 we naively expect to get 0.25, 
0.5 and 0.25 for the branching fractions to e~^e~ , e/j, and respectively, a consequence 

of the fact that the tau will go as often to a as an e. As can be seen in table 4 this is 
what we get at SPSlb while SPSla gives slightly larger values for the same-flavour (ee and 
HIJ.) channels. This difference is consistent with contamination by lepton pairs from the 
cascade decay chain, which is important for SPSla but not for SPSlb. 

Since this kind of contamination will be most important for tau-rich operators where 
there are fewer leptons from the neutralino decay, we also expect to see it in L1L3E3, which, 
in fact, we do. For a L1L3E3 coupling, we expect the same amount of e'^e^ and en from 
the er channel and from the rr channel we expect contributions to e+e", e/j, and n'^ 
in the ratios 1:2:1 (again due to the two possibilities for e/x). The result should be that 
the branching fraction to e/x is larger than the fraction to e"'"e~ by about the same amount 
as the total branching fraction to and this is rather close to what we see for SPSlb, 

but at SPSla the e^e" and e/x values are of the same size, indicating some contamination 
in the e"'"e~ (and jJL^ jJL^) channel from the cascade chain. 
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Four out of the nine LLE operators are such that they ahow one channel with prompt 
Ught leptons and one channel with a light lepton and a tau. These operators are LiL^Ei, 
L1L3E2, L2L3E1 and L2L3E2. Let us look at L2L3E1 as an example (the same logic 
applies to all of them). In this case, the prompt lepton channel is efx and, in addition, we 
get decays to er. Since the tau will decay equally often to electrons and muons, we expect 
a large signal in the channel both from the prompt leptons and the er channel. In 
addition we expect some contribution to the e~^e~ channel from the er decays, while /U"*"/!" 
should give nothing. As can be seen from tables 3 and 4 this is, to a good approximation, 
indeed the case. We also notice that the dominance of the largest channel is quite similar 
for all the SPS points considered and all the four couplings. However, there is a clear 
difference between the two ways of measuring the branching fractions; for L2L3E1 we have 
Pen ~ 0.81 while J Me^i ~ 0.87 and this difference is nearly independent of the choices for 
the SPS point and the R-violating coupling as well. 

The estimate of the branching fraction of the dominant channel will depend strongly 
on which cuts are used as well as how many leptons not originating from neutralino decay 
we have. Therefore it is very difficult to make any precise predictions, but the difference 
seen between the two measures is expected due to the differences in the cuts used. 

If we now compare the predictions from Eq. (3.8) with the direct measurement of the 
fraction (3.7), wc sec that operators with an er channel give a measured fraction about 
twice as large as the prediction (3.8) (see L1L3E1, LiL^E^ and L2L2E1 in table 3), while 
for operators with a fir channel the prediction (3.8) overestimates the fraction by about a 
factor two (see LiLsE2, L2LSE2 and L2L^E^ in table 3). This is in accordance with our 
expectation and illustrates that we can infer the presence of tau in the operator without 
having to rely on tau tagging. 

3.7 Redefining taus 

Another way of getting a grip on the decay channels containing taus is to include hadronic 
taus in the analysis; in this section leptons also include rjct's and the only cuts imposed on 
the leptons are the acceptance cuts of > 5 GcV (10 GeV for rjct) and isolation. 

The event selection used here is exactly four leptons (including rjet's), two positive and 
two negative. In addition, we require p^^^^ > 100 GeV. 

We can then expand the definitions (3.3) and (3.5) to include rjet's. Using these 
probabilities we can estimate all branching fractions of neutralinos decaying to lepton 
pairs. The result for all couplings at SPSla is shown in table 5 along with the result for 
L2L3E1 at SPSla, SPSlb and SPS6. It is clear from table 5 that this method does pick 
out the correct channels as the dominant ones; for all couplings, the two dominant ones are 
the same as those given in table 1. It should be noted that channels with tau will give rise 
to several different final states depending on the decay of the tau. If we look at L2L3E1 
at SPS6 the above points become clear; the two dominant final states are e/i and erjet and 
the er channel contributes about 0.05 — 0.10 to both Pg+e- and Pe^ and about 0.3 to -Perjet 
which is consistent with the branching ratios of taus. We also have some contamination of 
around 0.02 in all channels. 
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model 






p , 

'.lot '.ict 


p 


p 


p 


L1L2E1, 


SPSla 


0.38 


0.02 


0.004 


0.41 


0.08 


0.03 


L1L2E2, 


SPSla 


0.02 


0.40 


0.004 


0.40 


0.03 


0.07 


L1L2E3, 


SPSla 


0.08 


0.09 


0.06 


0.16 


0.30 


0.31 


LiL^Ei, 


SPSla 


U.44 


U.UOo 


O.Oo 


0.10 


0.36 


0.02 


L1L3E2, 


SPSla 


0.02 


0.09 


0.03 


0.44 


0.04 


0.33 


L1L3E3, 


SPSla 


0.09 


0.02 


0.26 


0.10 


0.40 


0.13 


L2L3E1, 


SPSla 


n no 

u.uy 


U.Uz 


n HQ 
U.Uo 


n A A 
U.44 


U.oZ 


U.Uo 


L2LSE2, 


SPSla 


0.005 


0.45 


0.03 


0.10 


0.02 


0.35 


L2L3E2, 


SPSla 


0.02 


0.09 


0.26 


0.10 


0.13 


0.39 


L2L3E1, 


SPSla 


0.09 


0.02 


0.03 


0.44 


0.32 


0.05 


L2L3E1, 


SPSlb 


0.09 


0.01 


0.04 


0.44 


0.35 


0.04 


L2L3E1, 


SPS6 


0.10 


0.02 


0.01 


0.54 


0.32 


0.02 



Table 5: Branching fractions for different final states, estimated from events with 4 leptons (in- 
cluding Tjet) with a total charge-sum of zero. 



Compared to the cases of table 3, there are now more leptons that do not arise from 
neutralmo decays. For example, for the SPSla point, each of the couplings considered 
results in more rjet's than expected from the neutralino decays alone. And when comparing 
the values for L2L^Ei at SPSla, SPSlb and SPS6 we notice some differences which are 
consistent with the differences in the contributions from the cascade decay chain; e.g. we 
see that SPS6 gives fewer Tjct's than SPSla and SPSlb. The conclusion is that the numbers 
are rather sensitive to the parameter point. 

These issues make it hard to interpret the exact values of the branching fractions. It 
is also important to remember that tau tagging efficiencies should be added on top of this 
and that will further suppress the rjet's in comparison with the light leptons. 

3.8 Combinations of couplings Case study 

To illustrate the usefulness of the observables discussed in the previous sections, let us 
study two cases where two operators are comparable, namely, 

mixl : L1L2E1 = L2LSE2 = 10-^ 

mix2 : L1L2E2 = L2L3E1 = 10"^. ^ ' ' 

These combinations are chosen such that the total flavour content of mixl and mix2 are 
the same at the operator level. However, the effect of the cuts will be somewhat different 
and that will induce some differences. For example, the pT spectra for the two mixtures at 
SPS6 (as shown in Fig. 7) do look similar, but we can see some small differences caused by 
the fact that the light lepton flavour accompanying the tau gets suppressed by the cuts. 

The corresponding invariant mass distributions (Fig. 8) are clearly different. For mixl 
we see prompt di-lepton (i.e. Eq. (B.2)) signals in each of the ee, e/i and /i/i channels. 
We also see that at masses below some 100 GeV, the latter two distributions are slightly 
larger than the ee one (we ignore the triangular contributions from the cascade chain). 
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Figure 7: Electron, muon and rjct spectra for mixl (left) and mix2 (right) at SPS6. 
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Figure 8: Di-lepton invariant masses for mixl (upper panels) and mix2 (lower panels) at SPS6. 

A fit confirms tliat the difference is well explained by an additional contribution of type 
(B.ll), which in turn suggests a /ir channel. It should be noted that the ability to detect 
a fiT channel in the light di-lepton distributions depends on the SUSY model so it cannot 
always be done, but we see that the //r channel clearly shows up in the //rjet distribution 
as well. 

Looking at mix2 in the lower panels of Fig. 8, we see a strong prompt e^u signal as well 
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as a smaller prompt fj,^ signal. This is difficult to explain without assuming at least two 
couplings where both give e/j, and only one gives n/j,. We also see a small signal in the ee 
distribution, consistent with Eq. (B.ll) and the fit also prefers such a contribution to the 
e/i channel. This suggests an er channel, which is confirmed by the erjet distribution. 



model 
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(3.7) 


(3.8) 


JMee 






mixl, SPSla 


0.27 


0.37 


0.36 


0.48 


0.84 


0.29 


0.35 


0.36 


mixl, SPSlb 


0.27 


0.37 


0.37 


0.44 


0.82 


0.29 


0.36 


0.36 


mixl, SPS6 


0.28 


0.37 


0.36 


0.48 


0.85 


0.30 


0.34 


0.37 


mix2, SPSla 


0.07 


0.62 


0.30 


1.31 


0.63 


0.05 


0.63 


0.31 


mix2, SPSlb 


0.08 


0.62 


0.29 


1.47 


0.65 


0.05 


0.64 


0.31 


mix2, SPS6 


0.08 


0.62 


0.31 


1.35 


0.63 


0.06 


0.62 


0.32 



Table 6: Similar to table 3 for mixl and mix2 at SPSla, SPSlb and SPS6. 
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mixl, SPSla 
mixl, SPS6 


0.19 
0.23 


0.24 
0.29 


0.02 
0.003 


0.25 
0.29 


0.05 
0.01 


0.21 
0.17 


mix2, SPSla 
mix2, SPS6 


0.06 
0.06 


0.21 
0.25 


0.02 
0.006 


0.42 
0.51 


0.19 
0.16 


0.06 
0.02 



Table 7: Similar to table 5 for mixl and mix2 at SPSla and SPS6. 



Let us now look at the counting observables from the last section. They are all listed 
in tables 6 and 7. If we first look at table 6, we again see that these numbers seem largely 
independent of the SPS point. When looking at mixl we see that we have significant 

branching fractions to ee, e/i and /i/i and that the last two are a little larger than the first. 
Wc also see that the measurement of the fraction (3.7) is about half the prediction from 
(3.8), suggesting a /xr channel which would also explain the higher values for e/i and /i/t 
compared to ee. 

If we instead look at mix2 we see a large branching fraction to e/t and also a signif- 
icant /t/( channel while the ee channel is small but non-zero. We also notice a measured 
fraction (3.7) that is twice that of the prediction (3.8), indicating an er channel, consistent 
with the small ee contribution. 

If we finally turn our attention to table 7, we see a larger variation with SPS point but 
the structure of the numbers is in agreement with our previous conclusions. 

To summarize, the invariant mass distributions and the counting observables from the 
last section can independently help identify the relevant decay channels which then leads 
us to the correct combinations of operators. The comparison between the measurement of 
the fraction (3.7) and the prediction (3.8) is the most reliable way of detecting an er or /tr 
channel without relying on tau tagging. 

4. LQD coupling 

LQD couplings can induce two types of neutralino decays. The neutralino can either go 
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to one up-type and one down-type (anti)quark and a charged lepton or to a down-type 
quark-antiquark pair and a neutrino. In the case of LiQ^D^ the latter decay dominates 
since otherwise a top would have to be produced, which is either strongly suppressed by 
phase space or requires an off-shell top; hence, in these cases, there are no isolated charged 
leptons from the decay and the lepton flavor of the A' operator cannot be determined. 

Table 8 summarizes possible scenarios for the various flavour combinations. Operators 
easier to study are of the type Li^2Qi,2Dz, where lepton identification combined with b- 
tagging might give some interesting signals. The reason the £j signal is mentioned in 
Table 8 is that in many cases the two jets can be so close in the detector that they merge 
into one jet and the peak expected in the ijj distribution instead shows up in ^j. This is 
also partially true for ih^et't however, in that case a more realistic algorithm for b-tagging 
is needed to know how well the b-tagging works if the jets are so close that they might 
merge. 

When the A' operator has r 
flavour the situation is rendered 
more complicated on account of the 
loss of information involving the 
neutrinos in the decay of the tau. 
The most problematic situation, how- 
ever, appears when we have a Q3 
operator, as discussed above we only 
expect decays to a neutrino and 
two jets and therefore this scenario 
might be hard to identify. 



Coupling 


Decay products 


Signals 


-^l,2<5l,2-Dl,2 


Iqq, uqq 




-^^1,2^1,2-03 


iqb, vqb 


ibjct, (^bjctj 


-^3Ql,2-Dl,2 


rqq, uqq 


Tied, Tjedj 


^3Ql,2-D3 


rqb, vqb 


Tjet^jet, Tjet^'jetJ 


^l,2,3Q3i^l,2 


vbq 




-^1,2,3Q3-D3 


ubb 





Table 8: The various possible scenarios with LQD opera- 
tors. It is here assumed that the neutralino is lighter than 
the top. Notation: rjet, &jet and j denote a tau jet, a b-jet 
and a light-flavor jet, respectively. 



4.1 Event selection 



The task of suppressing backgrounds is a bit more difficult with LQD couplings as compared 
to LLE couplings, especially when it comes to suppressing tt events. However, the lepton 
multiplicity is useful in this respect. 

To cover most possibilities we use three different cuts in this analysis. All three require: 

• Transverse sphericity > 0.2, 

• Et = ^leptons,jets\PT\ > 1000 GeV, 

• Hardest jet pr > 300 GeV. 



We define a "hard cut" to additionally require at least two same-sign isolated leptons 
with pt > 20 GeV. The same-sign requirement suppresses the ti background by several 
orders of magnitude, whereas the signal is reduced to only about 1% of its precut value. 

An alternative is to use a "loose cut" which passes a lot more of the signal, at the 
cost of a significant background from tt events. This requires only one isolated lepton with 
PT > 20 GeV, and at least 7 jets with pr > 20 GeV. 
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The third cut is meant for the cases where we only get neutrinos and no leptons. 
Therefore it requires instead of leptons, p^'**** > 200 GeV, and at least 7 jets with px > 20 
GeV. We name this the "pTmiss cut" . 



model 


hard cut 


loose cut 


pTmiss cut 


L1Q1D2, SPSla 


0.011 


0.082 


0.043 


L1Q1D2, SPSlb 


0.043 


0.24 


0.16 


LiQiL'2, SPSS 


0.06 


0.25 


0.16 


L1Q1D2, SPSS 


0.015 


0.089 


0.046 


L1Q1D2, SPS6 


0.008 


0.096 


0.083 


L1Q1D2, SPSla 


0.011 


0.082 


0.043 


L1Q3D3, SPSla 


0.0013 


0.03 


0.048 


L3Q1D3, SPSla 


0.0027 


0.043 


0.041 


tt 


1.5 • 10-6 


0.002 


0.00083 



Table 9: The fraction of events that passes the various cuts for some LQD operators, as well as tt 
production. Note that this does not take the different cross-sections into account. 

The effects of these cuts, on some representative LQD couplings as well as on the tt 
background, are summarized in Table 9. All other backgrounds are less important than ti 
for these cuts. 

We see that for the hard cut we get a satisfactory suppression of the ti while the loose 
cut will suffer from a significant contamination by tt events. 

It is also clear that the operators with tau flavour are more problematic due to fewer 
leptons. As expected, of the operators in Table 9, the most problematic one is LiQ-^D-^ 
where the signal would be i^ebb or e^tb. With the t-quark above the kinematical threshold, 
only the t'e66-channel is open and we have to rely on the pTmiss cut. This last case will 
appear very similar to an R-parity conserving scenario, the most significant difference being 
that the particles now responsible for the missing pr are the (massless) neutrinos rather 
than the neutralinos. This problem is present for all operators of type LQ^D. 

4.2 Invariant mass distributions 

There are plenty of interesting invariant mass distributions to be studied with respect to 
LQD operators. Unfortunately, the lack of knowledge of the charge of light quark jets 
implies that the same-sign subtraction tool is essentially rendered inoperative. As a result, 
most signals tend to drown in combinatorial and other backgrounds. This problem would 
be reduced if the couplings were small enough for displaced vertices to appear, however, 
as mentioned before, we do not make use of such possibilities here. 

A further problem is that the boost of the neutralino results in all the decay products 
being very close in the detector, thereby raising the probability of jets merging. It is, 
therefore, not always true that we get two jets and one lepton from each neutralino. On 
the other hand, the closeness of the decay products can be used to choose which jets and 
leptons to compute the invariant mass of. 
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Figure 9: Invariant mass distribution for L1Q1D2 at SPSla: upper left: violet: £j distribution, 
black: ijj where the jets are closer than AR = 1 to the lepton, this distribution has been multiplied 
by 5 to be visible. Upper right: ijj distributions using Rj^t = 0.2. Lower left panel: same as upper 
right but with _Rjct = 0.4, lower right: same as lower left but for SPSlb. In this figure, the hard 
cut has been used. 



These proximity issues make the choice of jet algorithm very important. Especially 
the parameter i?jet is very important. To optimize each analysis, we perform it using jets 
as defined by Rjet equal to 0.2, 0.4 and 0.7 and show results using the most favourable 
value for each analysis. The event selection always uses jets as defined when i?jet = 0.4. 

In Fig. 9 these points can be seen for L1Q1D2 at SPSla: if ijj where the two jets are 
within AR = 1 from the lepton, are used, one can get a very nice peak at the neutralino 
mass (black curve, upper left panel). If all ijj are used, the peak is less pronounced but 
can be clearly seen both with i?jet = 0.4 as in the lower left panel, and even better when 
Rjet = 0.2, which can be seen in the upper right panel. Since the two quark jets quite 
often merge into one, it is also possible to see a nice peak at the neutralino mass in the 
distribution of all ij combinations, as seen in the upper left panel. It is also worth noting 
that the branching fraction for neutralino decay to ijj as compared to the fraction to z^jj 
depends on the SUSY parameter point and e.g. SPSlb has a more favourable value, which 
can be seen in the lower right panel of Fig. 9, where the peak is clearly more pronounced 
than in the comparable plot for SPSla (lower left panel). 
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So far wc have discussed the case with one hght lepton and two light jets in the decay. 
Another group of operators have b-quarks in the decay and that opens up for b-tagging. 
Assuming that the b-tagging can also determine the charge of the 6jet , it becomes possible 
to use same-sign subtraction to reduce the combinatorial background. 

Fig. 10 shows Ibjetj distributions as well as ibjet distributions where same-sign sub- 
traction has been employed on the £bjet pair. In the upper left panel we see that this gives 
a good signal for LiQiAs where we again have a peak at the neutralino mass. The signal 
in the e6jet channel looks like the distribution (B.2) with a small peak at the high end 
(caused by merged jets), which is consistent with our expectation. Due to the low number 
of charged leptons in the neutralino decay at SPSla, we have used the loose cut to achieve 
the upper panels of Fig. 10. For SPSlb we get more charged leptons and as is clearly seen 
in the middle panels the hard cut is fine for extracting these signals. However, the SUSY 
production at SPSlb is also about one order of magnitude smaller than at SPSla so using 
the hard cut is more important to reduce the backgrounds. 

If the lepton is a rjet, things are much more difficult; first of all we lose energy to the 
neutrino, so we do not expect a peak anymore. However, if we use Rjet = 0.2 together with 
the pTmiss cut for SPSS we see that the expected signals are clearly there. In the rjetftjeti 
distribution (lower left panel) we see a bump below the neutralino mass. Perhaps the best 
and most conclusive signal though, is in the Tjctfcjct channel (lower right channel) where 
we see a distribution similar to what we discussed in appendix B for tau decays (the ir^et 
distribution of Fig. 21). 

4.3 Mixing LQD operators 

The chances to identify two or more simultaneously large LQD type operators, depend 
strongly on the specific combination; if the two operators are identified through different 
channels, the identification can be rather straight-forward. 

One example of this can be seen in Fig. 11 where one operator gives a signal in /ubjcti 
and //6jet channels and the other shows up in the ejj distribution. This case is rather 
fortunate since the two operators can be identified independently; if the lepton flavour of 
the operators would have been the same, it would have been much harder to tell that we 
have one operator with and one without a component. 

In Fig. 12 we show another combination where an operator (A211) giving a peak in 
the njj distribution, is combined with a more problematic tau-fiavour coupling (A3^3) and 
again we see that both operators can be identified successfully. 

In conclusion, having more than one large LQD operator at the same time is no 
problem as long as the signals from the couplings appear in different channels and the 
operators are identifiable on their own. 

4.4 Backgrounds 

As mentioned before, the most important background is tt. Although the hard cut is very 
effective in suppressing this background; the need to sometimes use the loose or pTmiss cut 
makes it necessary to look at the expected invariant mass distributions from the tt events. 



-25- 



0.015 



0.01 



0.005 



■SPSIa 




— eb] 

— Vb] 






— xbj 




J 




.... 1 . 


X 




1 ........ . 



50 100 150 200 

M[GeV] 



0.004 



0.003 




50 100 
M[GeV] 



150 



200 



0.015 



0.01 



0.005 




0.004 



50 100 150 

M[GeV] 



200 




50 100 
M[GeV] 



150 



200 



0.015 



0.01 



0.005 - 




50 100 
M[GeV] 



200 



0.004 



0.003 



0.002 




50 100 
M[GeV] 



150 



200 



Figure 10: Left panels: ibjctj distributions with same-sign subtraction in the ibjct pair. Right 
panels: ibj^t distributions with same-sign subtraction. Upper and middle panels: LiQiD^ at SPSla 
(upper) using the loose cut and SPSlb (middle) using the hard cut. Lower panels: L3Q1D3 at SPSS 
with i?jct — 0.2 and using the pTmiss cut. 



The interesting distributions here are the invariant masses including 6jet since that is 
where the weaker signals occur and we expect some contribution there from tt due to tops 
decaying to ibjeti^)- 

These distributions are shown in Fig. 13 where we can see that the only significant 
distribution is in the fibjet and eftjet channels for the loose cut. The distribution can be 
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Figure 11: Combinations of LQD operators: Invariant masses for \'^2i = 10^'' and A223 ~ li)^^ 
at SPSlb. The two operators differ both in lepton flavours and in quark flavours. Upper left: Ij 
and tjj (magnified by a factor 5), in the latter case the two jets are close to the lepton (see Fig. 9). 
Upper right: Ij]. Lower panels show Ib-^ctj (left) and ib-^^t (right) distributions. All plots use the 
hard cut. 



understood as a triangular distribution ending around 150 GeV from the decay of the top 
and a combinatorial background that due to the same-sign subtraction and the charge 
correlation between the top and antitop, is negative. 

It is also important to note that we see no significant signal when the pTmiss cut and 
i?jet = 0.2 are used. Simulations where standard model processes were added, also show 
that the studied distributions are not significantly changed by these backgrounds. 



5. UDD coupling 

Of the various scenarios with trilinear R-parity violation, the one that is most difficult to 
analyse at a hadronic collider is that of non-zero JJDD couplings. With the lightest neu- 
tralino now decaying into three jets, extracting the signal from the large QCD background 
is very challenging [37]. In fact, even for simpler situations like resonance production or 
one-step decays, the reach of the Tevatron or the LHC is very limited [15,38]. 

As mentioned earlier, assuming moderately light neutralinos, we expect them to be 
highly boosted due to the large center of mass energy of the collision. This means that the 
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Figure 12: Same as Fig. 11 for = 10 ^ and Agj^g = 10 ^ at SPSlb. Here the two lower panels 
use the pTmiss cut as well as i?jet = 0.2 to get enough Tjet's. 



decay products of the neutralino will be rather collimated, which has enormous implications 
for the identification of LQD operators as discussed in the previous section. This would 
also be the case for UDD operators where the three jets from the decaying neutralino 
will sometimes be so collimated that they merge. This very feature can be turned to an 
advantage; a UDD signal may be extracted by studying jet substructure, where one looks 
for a big jet that looks like it could consist of three jets that have merged [39]. 

There is, however, an exception to the above picture, namely any of the three operators 
U^DiDj, for the top flavour of these couplings ensures a very different behavior. If the 
neutralino is lighter than the top, it may appear stable on the detector scale and we might 
be facing a fake MSSM scenario. On the other hand, if the neutralino is heavier than the 
top, it decays into a top (or antitop) and two jets that often would tend to be soft (due to 
limited phase space). While a tt + n-jets final state has a very large QCD background, note 
that the latter decreases fast with both increasing n as well as with an increasing p™™(jet) 
(the minimal requirement on the transverse momenta of these additional putative jets). 
And as additional jets accrue not only from the decay of the neutralinos, but also in earlier 
steps in the SUSY cascade, and as at least some of them are likely to have a sufficiently large 
Pt, this final state is worth investigating. However, rather than performing this admittedly 
difficult task, we recommend concentrating on a relatively cleaner subsample. Noting that 
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Figure 13: Invariant mass distributions for events. The left panels show ^6jotj distributions and 
the right panels show £bjet distributions. In the upper panels the loose cut is used and in the lower 
panels the pTmiss cut is used and i?jct = 0.2. 



the neutralinos are Majorana particles, the very last pair of decays is as likely to produce a 
like-sign top-pair {tt or tt) as an opposite-sign one. The possibility of same-sign tops makes 
the extraction of these events rather straight forward. One way is to look for same-sign 
isolated leptons accompanied by two same-sign b-jets; this approach was taken in a similar 
study [40]. 

However, it turns out that the harder cuts used on LQD operators are also efficient in 
this case. This comes as no surprise since those cuts are based on same-sign leptons. If we, 
for example, look at U3D1D2 at SPS6 where the neutralino mass is 189 GeV so that the 
top decay channel is open, about 1% of the events pass the hard cut from the LQD section. 
This leaves good hope for detection and if we then look at the ibjct invariant mass, as is 
done in Fig. 14, we see the same structure (although rather weak) as in the corresponding 
distribution from tt events as shown in the upper right panel of Fig. 13. This similarity 
demonstrates the presence of top quarks in the events and can therefore be used to confirm 
the presence of a U^DjDf^ operator. 
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Figure 14: Invariant mass distributions for \'^i2 SPS6. The distributions shown are £6jot 
distributions for ah lepton flavours. 

6. Mixing and models 

To date, most phenomenological analyses have assumed the dominance of a single R- 
violating operator. However, once a flavour structure is invoked (e.g., to explain fermion 
masses and mixings) [41], hierarchies in R- violating couplings would be related to the 
flavour charges of the respective fields [21]. Even more importantly, even if only one R- 
violating operator is to be postulated in the interaction basis, fermion mixing would, in 
general, induce non-zero values for others. This implies that it is natural to expect a range 
of hierarchies for R-violating operators. 

6.1 LLE and LQD operators 

Before entering into a discussion of models that require operators being related in mag- 
nitude, we shall first look at some examples of the presence of two different operators. 
Without jeopardizing proton stability, we can allow for the simultaneous presence of LLE 
and LQD operators. 

Since both LQD and LLE operators can be identified due to the lepton multiplicity, 
the presence of one might mask the presence of the other, especially if the corresponding 
decay widths are comparable. Interestingly, for many of the SPS points in question, having 
one each of LLE and LQD operators of the same magnitude, results in the decay width 
associated with the LLE coupling being about 10 — 20 times larger than that due to the 
LQD operator.^ Furthermore, even if the decay widths for the two sets of channels were 
comparable, the LLE ones would be more visible due to their higher lepton multiplicity 
which implies easier detection. Thus, in order to study scenarios where one might actually 

*The reason for this is not far to seek. The squarks tend to be much heavier than the sleptons, and, 
of the latter, the stau tends to be the lightest. Consequently, the purely leptonic decay modes of the LSP 
dominate over the semileptonic ones. 
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Figure 15: Mix of operators: Invariant mass distributions for A121 = 10^^ and A']^23 — 10^^ at 
SPSla. Upper panels: di-lepton invariant mass distributions. Lower left panel shows Ib-^etj and 
lower right shows £6jet- 



see effects of both operators at the same time, we take the LQD operator to be 10 times 
larger than the LLE one. 

While most of the conclusions from the study of single couplings are still valid when 
more than one coupling is relevant, the details of section 3.6 are modified. In particular, 
the differences between the measured and estimated values of the fraction (3.7) no longer 
point to the relevant tau channels. Therefore, it is more useful to study invariant mass 
distributions. 

A favourable scenario would be the combination of two operators that are both rela- 
tively easy to identify. One example of this would be A121 = 10~^ and A']^23 = 10"'^. The 
invariant mass distributions corresponding to this scenario, at SPSla, are shown in Fig. 15. 
We clearly see the signals expected from the LLE coupling in the and ee distributions 
(upper left panel) as well as the ^6jeti (lower left panel) signals expected from the LQD 
coupling. We also note that the distributions including rjet's (in the upper right plot) show 
signals that might suggest a rr channel, however, those tau pairs actually come from the 
cascade chain and one has to be careful not to interpret them as decay products of the 
neutralino. 

Whether the potentially confusing background from the cascade is present or not. 
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depends on the parameter point, in e.g. SPS6 it is not. However, SPS6 comes with another 
problem; it turns out that the branching fractions for decay to charged Icptons plus quarks 
is much smaller (factor Ri 7) than the branching fraction to neutrino plus quarks. This, 
then, results in fewer events with leptons and, as a consequence, the identification of the 
LQD operator is rendered more difficult. 

Similarly, for other experimentally complicated scenarios, such as operators with a lot 
of tau flavour, the problems (in identifying the LQD operator) noted in section 4 remain 
pretty much unchanged when combinations are considered. In general, the identification of 
the LLE and the LQD operators can be handled independently and the only interference 
between them is that for the smaller coupling we get less statistics to work with. This 
also means that there is an interval where both operators can be identified while in the 
periphery we will see only one dominant operator. This interval is centered roughly where 
the LQD operator is about 10 times larger than the LLE coupling. In general, the LLE 
operators are easier to identify, which allows for smaller LLE couplings to be successfully 
identified. 



6.2 Link with flavour symmetries: Two examples 

For a given flavour model, the relative magnitude of R-violating operators is, to a large 
extent, predicted, allowing us to tabulate the expected signals in neutralino decays. This 
will be done for two representative cases, namely Left-Right symmetric models and SU{5). 



6.2.1 Left-Right-symmetric models 



The simplest starting point is a single U{1) family symmetry with the same charges for the 
left- and right-handed states (left-right symmetry) as shown in Table 10, where, e.g., the 
choice a,; = (—4, 1,0) [2] gives an acceptable pattern for the mass matrices. 

With these charge assignments the 
quark mass matrices (up to numerical 
factors and phases, which, in general, 
are expected to be of order unity) take 
the form 



C/(l) 



Qi Ui Di Li Ei H2 



Hi 



hi -2a3 -2a3 



Table 10: Assignments ofU{l) charges. 



/e8 e3 e 



e3 e2 e 
e 1 



M' 



down 



e 1 



(6.1) 



where e y/e ^ 0.2. 

We now consider the effect of the U (1) symmetry on the pattern of allowed -R-violating 
interactions [21]. In this simple example, with all fermions of a given family having the 
same charge and with a left-right symmetry, the charges of the operators depend only on 
the combination {i,j,k) and are independent of the type, viz. LLE, LQD or UDD, (see 
Table 11). 

The parameter w accounts for the fact that the charge assignment is not unique in 
model constructions [22] (although it is strongly constrained by phenomenological and 
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ijk 


111 


121 


122 


222 131 


m 


-12 - w 


— 7 — w 




3 — w —8 — w 


ijk 


133 


333 


223 


233 123 


U{1) 


-4-w 


—w 


2-w 


1 — w —3 — w 



Table 11: Operator charges in a model (see text) with both family and Left-Right symmetry. Here 
w parametrises flavour-independent contributions [22]. 

theoretical arguments). For instance, while the addition of flavour independent contribu- 
tions would not modify the fermion mass and mixing ratios, it can affect other consider- 
ations, such as anomaly cancellation conditions. Flavour-independent contributions could 
also arise from additional fields with a non-trivial flavour charge that couple to all opera- 
tors. These issues in the context of R-violating hierarchies are addressed in more detail in 
Ref. [22]. 

The above flavour symmetries cannot ensure, by themselves, that rapid proton decay 
is avoided. This is done by imposing a baryon or a lepton parity [2,3,18], which imply that 
the cases with AL ^ and with /S.B ^ have to be considered separately. For instance, 
lepton- number violating operators can be eliminated by imposing a lepton triality [18], 
under which the flelds transform as 

Z3 : (Q,C7,L>,L,^,ifi,J/2)^ (l,l,l,a,a',l,l). (6.2) 

This allows only the baryon-number-violating operators and the mass terms, while forbid- 
ding lepton-number-violating ones. On the other hand, to forbid baryon-number violating 
operators, we would work instead with a baryon triality, such as in Ref. [2], viz. 

Z3 : {Q,U,D,L,E,Hi,H2) (1, a^, a, a^, a^, a^ a). (6.3) 

Along these lines, we can use the hierarchies implied by Tabic 11 to make the link 
between flavour models and the expected decay rates in neutralino decays for two separate 
cases, namely: 

• for LLE and LQD operators; 

• for UDD operators. 

While w can be adjusted in order to ensure that all operators remain within the experi- 
mental bounds, there can be additional sources of suppression in the couplings (arising due 
to a small tan P in supersymmetric models, the form of the Kahler potential, or additional, 
model dependent, features of the theory that may involve extra fields and symmetries). To 
keep the problem manageable (especially since our purpose here is to show how specific 
models can be tested in fiavour neutralino decays), we will: 

• assume that the couplings that correspond to higher fiavours are bigger (have smaller 
flavour charge). The Yukawa couplings generating fermion masses are larger for the 
higher generations, having a smaller net flavour charge. If this is true for R-violating 
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couplings as well, then the operators that involve third-generation flavours would 
dominate also in this case. Dominance of the R-violating couplings of the heavier 
flavours would also favour radiative over three body gravitino decays in scenarios 
with gravitino dark matter [17,18]. 

• ensure that all couplings remain within the experimental bounds. 

For AL 7^ the strictest bounds are on LiQiDi from nuclear f3f3 decay and on LiL^E^ 
from bounds on Majorana neutrino masses, constraining the choice of the charge w to be 
|12+-u;| > 2 and 14+^)1 > 2, which are easy to satisfy [22]. We also note that the magnitudes 
of the couplings in Table 11 are symmetric in the three indices ijk. This implies, for 
example, that, at this level, the X'121 and X'112 couplings should have similar magnitudes. 
This must be made consistent with the constraint {L1Q2D1) • {L1Q1D2) < 4 • 10~^, which 
arises from bounds on Attik- In the present context, this constraint indicates that the 
relevant charge |7 + has to be large. 

The problem gets more complicated by taking into account mixing effects, which for 
left-right symmetric models have been studied in detail in [22] , in combination with bounds 
on both individual couplings as well as on products. Summarising, it turned out that the 
strong correlations were leading to a suppression of all couplings, and even the 333 flavour, 
in order to be fully safe, had to be smaller than 0.006. Within this framework, therefore, 
single superparticle productions are suppressed and the best signal would be pair pro- 
ductions followed by i?-violating decays (which favour the neutralino search channel even 
further). Finally, quantum corrections, expressed through renormalisation group effects, 
would only mildly affect the relative hierarchies of R-violating couplings (although they 
would modify their absolute values between the GUT and the low energy scale) [9]. 

For the AL 7^ operators, we may then study the correlated operators: 

LLE^s + LQDeff, (6.4) 

with 

LLE^s = (ML-^E^ + e^LiL2E2 + 0(e'^), 

LQD^s = L3Q3D3 + e{L2Q3D3 + LgQs^'a + ^^3(33^2) + 0{e^), (6.5) 

Due to the larger decay widths associated with the LLE operators as compared to the 
LQD operators as well as the lack of a clear signal associated with the Q3 factor in the 
dominant operator; the only thing wc can sec is the L2L3E3 operator and possibly a hint 
of the L1L2E2 operator. This is illustrated in Fig. 16 where the leptonic invariant mass 
distributions are plotted. The data is consistent with a dominant contribution from a /xt 
decay channel giving the //rjet distribution, together with contributions of type (B.ll) to 
the fijj, and e/x distributions. Since the fx/j, and e/j, distributions have kinematical cutoffs of 
the same scale as the /urjet distribution, we need an additional contribution of type (B.2) 
to those distributions, consistent with a subdominant L1L2E2 operator. 

Due to the dominant LQD operator, a larger fraction of the leptons in the final state 
will come from the upper decay chain rather than the neutralino decay, as compared to the 
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Figure 16: The same as Fig. 5 for the Left-Right symmetric scenario defined by (6.4) at SPSIa. 

pure LLE case, and therefore the numerical observables of section 3.6 are too dependent 
on SPS point and not very useful. This contamination also causes the invariant mass 
distributions to be much less smooth and therefore more difficult to interpret. 

If we now look at the UDD operators, we see that even though, within this model, 
the largest coupling is U3D2D3, the phase-space suppression caused by the top quark in 
the final state makes the dominant channels the ones with three quark jets, even at SPS6 
where the neutralino is sufficiently heavy to decay to a top. 

6.2.2 5?7(5) 

Another interesting possibility is that the family symmetry commutes with an SU (5) GUT, 
where the SM fermions are assigned as follows to the representations of the group: 

Q(i,d^). = Ql (6.6) 

From the above it immediately follows that : 

1. The up-quark mass matrix is symmetric (both left- and right-handed up quarks are 
in the 10, and thus have the same flavour charge). 

2. the charged lepton mass matrix is the transpose of the down quark mass matrix. 

In this case, a viable choice of charges obeying the restrictions of the symmetry (see e.g., 
[42]), is: 

Ql,2,3 = -£'1,2,3 = 3,2,0 

^1,2,3 = ^1,2,3 = 1,0,0 (6.7) 

leading to matrices that, apart from other features, lead to a maximal 2-3 lepton mixing. 
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where e « 0.2. 

Let us first look at the imphcations for the LLE operators. Since the charges of L2,3 
are the same, couplings such as LiL2Ek and LiL^Ek would be expected to be of similar 
magnitude. In short, the U{\) assignments of Eq. (6.7) lead to operator charges as listed 
in Table 12. 



ijk 


121,131 


231 


122,132 


232 


123,133 


233 


U{1) 


A — w 


3 — u; 


3 — u; 


2-w 


1 — w 


—w 



Table 12: LLE charges in SU(5) enhanced by a U(l) flavour symmetry. 



Similarly for LQD, where now we have the connection A^^-j, = Ajj^ arising directly from 
the way wc accommodate the fields in the GUT representation. For the LQD operators 
we also note the following: 

• Since the ^7(1) charges of ^2,3 are the same, the respective operators are linked, 
reducing the number of independent couplings. 

• Since the U{1) charges of ^2,3 are the same, the respective operators are linked, 
reducing the number of independent couplings. 



This leaves us with the results of table 13. Finally, for the baryon-number violating oper- 
ators, we find the results given in table 14. 



ijk' 


111 


112,113 


121 


122,123 


131 


132,133 


Uil) 


5 — w 


A — w 


A — w 


3 — w 


2-w 


1 — w 


ijk' 


211,311 


212,213,312 


221.321 


222,223,322,323 


231 


232,233.332.333 


m 


A-w 


3 — w 


3 — w 


2-w 


1 — w 


—w 



Table 13: LQD charges in SU(5) enhanced by a U(l) flavour symmetry. 



ijk" 


112,113 


123 


212,213 


223 


312,313 


323 


U{\) 


A — w 


3 — w 


3 — w 


2-w 


\ — w 


—w 



Table 14: UDD charges in SU(5) enhanced by a U(l) flavour symmetry. 



As for the L-violating part, for = 0, the L2L3E3 coupling is as large as the largest 
of the LQD couplings and therefore the prospects of accurately identifying this coupling 
are improved. Moreover, the leading subdominant LLE operators, namely LiL^E^ and 
L1L2E3, also carry a lot of tau flavour. This ensures an even larger dominance of the 
L2L3E3 operator after cuts are applied (in the Left-Right symmetric case the cuts partially 
compensate for the difference in width between the L2LSE3 and the L1L2E2 operator). All 
in all, the prospects for identifying the L2L3E3 operator are good but it does not appear 
possible to extract any of the subdominant couplings. 

This can be seen in the invariant mass distributions shown in Fig. 17 where all distri- 
butions look much like we expect for a pure L2L3E3 operator. Note that all the light lepton 
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distributions of the left panel of Fig. 17 vanish between « 80 GeV and the kinematical 
end-point of the //rjet distribution ^xl^ therefore exclude any significant 

neutralino decay to light leptons. 

Note that this strong statement is linked to the simplifying assumption that like for 
the Yukawa couplings, R- violating couplings of the third generation are in general larger. If 
this assumption is dropped, the picture would change. For instance, for the LLE operators: 

• w = 1 would imply a simultaneous dominance of the L1L2E3 and couplings, 

• w = 2 leads to a dominant L2L3E2 coupling, and 

• w = 3 results in a simultaneous dominance of the L2L3E1, L1L2E2 and L1L3E2 
couplings. 

A similar situation (but for different respective flavours) would occur for the rest of the 
operators. 

For the baryon- violating case, w = leads to a stronger dominance of the U^DjDk 
operators as compared to the Left-Right symmetric case. Consequently, the decay to top 
plus 3 jets will dominate if the neutralino is heavy enough, e.g. in SPS6. 

7. Summary 

Within a super symmetric theory with trilinear R-parity violation and a neutralino LSP,^ 
one would expect the usual pair production of squarks and gluinos, followed by cascade 
decays down to the neutralino, to be the dominant scenario at the LHC. The novel feature 
introduced by the R-parity violation would then be the three-body decay of the neutralino. 

Within such a model, the neutralino can decay through any of the 45 trilinear R- 
parity violating operators and therefore this channel allows us to study all those operators 
and their internal hierarchies simultaneously. On account of the typically large number of 

^The lightest neutralino could even be the NLSP if, for example, we invoke a gravitino LSP for dark 
matter; such an assumption does not affect this study. 
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detectable particles (especially leptons) produced in such ncutralino decays, the inferences 
drawn from such a study are not as dependent on the underlying SUSY scenario as R-parity 
conserving SUSY scenarios can be. 

In general, bilinear and trilinear terms are generated in different ways (for instance 
an underlying string theory would favour trilinear couplings; on the other hand, flavour 
symmetries could in principle favour either possibility). See, however, Ref. [45]. The 
decay modes Xi ~^ W^f^ and x? ~^ Zu exist in bilinear schemes, leading to interesting 
signatures [46,47]. For these terms to dominate, one would require non-negligible neutrino- 
neutralino mixing. 

We have shown that the prospects of identifying trilinear operators of the LLE type are 
rather good, mostly due to the abundance of invariant mass distributions one can explore. 
While one can, in principle, also identify such operators from pure counting measures, the 
corresponding conclusions would not only be more sensitive to backgrounds, but also to 
potential competing operators, as well as to the SUSY model. In addition to a significant 
reduction in such dependences, the invariant mass distributions are also promising for 
measuring the neutralino mass; this has been shown to be feasible thanks to the knowledge 
of the theoretically expected distributions, despite the lack of clear peaks or edges. 

As compared to LLE operators, LQD operators have both advantages and disadvan- 
tages. One disadvantage is that background in the form of tt events is more challenging, 
but the large lepton multiplicity is still sufficient to get a clean event sample. Operators 
that lead to charged leptons and jets also have the advantage of producing peaks in the 
relevant invariant mass distributions, which not only allow good determination of the op- 
erator, but also accurate measurements of the neutralino mass. Particularly promising are 
the operators leading to 6-quarks in the decay since, in this case, 6-tagging can also be used 
to extract the signal. 

With a tau as the lepton of the operator, things are not equally promising, due to 
the loss of energy to neutrinos in the tau decay. However, if the tau is accompanied by a 
6-quark, then b and tau tagging together might provide the information we need. For the 
case with a tau and two light jets, and for LQ^D operators (which only gives neutrino and 
two jets) we do not have an obvious method of identifying the operator. 

If more than one lepton number violating operators are large at the same time, we 
have shown that the prospects of identifying them are good (in most cases, the signals we 
are looking for appear in different channels for different operators and these signals are 
mostly independent of each other). It is also worth mentioning that if both LLE and LQD 
operators are large at the same time we expect to see only the LLE operator since it gets 
a larger branching fraction and is easier to detect. 

From flavour models one may expect large R-violating hierarchies, similar to those of 
the Yukawa couplings that generate fermion masses. Fermion mass terms are dominated 
by the heavier flavours. If this feature persists for the R-violating couplings as well, it 
would make the identification of the LQD operators very difficult since the dominant ones 
will either contain tau flavour or a operator, with either case implying a low chance 
of identification. In addition, accompanying LLE operators have the advantage of larger 
branching fraction as well as easier detection. The conclusion is that in this subclass of 
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Figure 18: Kinematics of three-body neutralino decay. 

flavour models we would expect to see dominant LLE operators with heavy flavours, i.e. 
lots of tau flavour. 

Finally we have looked at the case where a heavy neutralino (M^o > Mtop) decays 
via a U^DjDis operator to a top quark and two jets, and have shown that this scenario 
can be identifled. This is important for flavour physics since in many models the R-parity 
violating couplings with heavy flavours will be the dominant ones. 
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Appendix A. Energy spectra from three-body decays 

For any decay into a flnal state comprising more than two particles, the energy distribution 
of a given flnal state particle is determined not only by kinematics, but has a dependence 
on the dynamics as well. In particular, let us consider the decay of a neutralino to three 
fermions, namely 

Xi ^ /1/2/3. (A.l) 

There are three contributions, characterized by the identity of the sfermion (correspond- 
ing to the fermion coupling to the neutralino). Adopting the notation of Fig. 18, the 
corresponding matrix element will be of the form 

M ~ [u{pi)u{p)] [u{p2)v{p3)]- (A.2) 

Here, we disregard the energy dependence of the propagator (valid for a heavy sfermion) 
and the pseudoscalar part of the coupling, since its presence will not modify any differential 
decay rate. Then, 

Y,^^^^(p-Pi)iP2-P3) (A.3) 

spins 

In the decay, two of the flnal-state particles will be charged, and potentially observable. 
For a particular flnal-state particle, the energy distribution depends on whether it comes 
from the upper or lower vertex in Fig. 18. With the notations p = {M, 0}, pi = {Ei,Pj^}, 
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and ei = Ei/E^'^ = 2Ei/M, we find 

ldr^(l2eKl-ei) i = l, 
Tdei \2ef(3-2ei) z = 2,3. 

For i = 1, the distribution vanishes at the upper end due to angular momentum conser- 
vation. For comparison, if the distribution had been determined by kinematics alone, it 
would have read 

To see the implications of this, consider the operator L1L2-E3. The last factor represents 
a coupling to a tau or a stau (both charged), whereas the first two (by SU(2) invariance) 
represent couplings to a charged field and a neutral one, one of which will be of the first 
family, whereas the other will be of the second family. We thus have two distinct decay 
channels: 

X? ^ r+e-i/;,, xi r+n-ue, (A.6) 

and charge conjugates. Each of these channels is governed by three diagrams of the type 
shown in Fig. 18. Assuming that all the sfermions have comparable masses, and neglecting 
interferences, the energy distribution for each of the decay products would be averaging 
the distributions (A.4) over all three values of i, resulting in 

IdV 2 20 



The above assumptions are not entirely justified, e.g. we do expect differences in sparticlc 
masses to introduce differences between the three diagrams due to propagator effects and 
we have no reason to ignore interference altogether. However, Eq. (A. 7) should still be a 
good approximation and the assumptions involved provide a useful formula. 
The distributions (A.4), (A.7) as well as (A.5) are shown in Fig. 19. 

Appendix B. Theoretical invariant mass distributions 

In order to calculate the theoretically expected invariant mass distributions [43, 44] from 
the leptonic neutralino decays, we start with the energy distribution of the decay products, 
as given by Eq. (A.7). 

Let us now denote the three decay products a, b and c, with corresponding four- 
momenta Pa, Pb and Pc- The invariant mass M^b of particles a and b is given by 

M!b = {Pa + Phf = m2 - 2ME, + ml (B.l) 

where Ec denotes the energy in the neutralino rest frame and M is the neutralino mass. 
Assuming that particle c is massless and follows the distribution given by Eq. (A.7) we get: 

fM^,{Mab) = ^iM'- Mlf{M^ + hMl,). (B.2) 
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Figure 19: Energy distributions for the decay products in a three body decay, Eq. (A. 4). The 
straight black line shows the distribution achieved when only phase space (PS) is taken into account. 
The average, Eq. (A.7), is denoted 'RPV. 
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Figure 20: Theoretical invariant di-lepton mass distributions for massless particles from Xi decay. 
Left: No cuts. The £i curve is given by Eq. (B.2). Right: Cuts as indicated by Eq. (B.ll). Solid 
curves: directly produced leptons; long-dashed curves: one lepton comes from a decaying tau; 
short-dashed curves: both leptons come from decaying taus. Blue heavy curves do not take tau 
polarization into account while the thin green curves do. The labels refer to the heavy curves 
(blue). 

This distribution is shown in Fig. 20 (solid curves). 

If one of the charged particles is a r, we get a contribution to the di-lepton invariant 
mass only if the r decays leptonically. To calculate the distribution from this, we need to 
estimate a distribution for the fraction of the energy of the tau that is transferred to the 
charged lepton. 

To do this, we assume that the tau decays into three massless particles (one charged 
lepton and two neutrinos). The kinematics of a tau decay will be the same as discussed in 
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Appendix A for a neutralino decay. However, with a V — A coupling, i = 1 in Eq. (A. 4) 
represents the antiparticle among the decay products. In the rest frame of the tau, the 
charged lepton will then follow the distribution for ? = 2, 3 with M = Mt- 

The next step is to perform a Lorentz boost to the rest frame of the neutralino. If 
we denote the momentum four-vector of the charged lepton in the rest frame of the tau, 
(El, Pi) and that of the tau in the rest frame of the neutralino, {E,P), we get the energy 
of the lepton in the rest frame of the neutralino as: 

EiE Pi P , 
— \ . (B.3) 

Mr Mr ^ ' 

Since we assume that the lepton is massless, i.e., Ei = \Pi\, and that the lepton carries a 
fraction x of Mr/2 in the rest frame of the tau, i.e. Ei = xMr/2, Eq. (B.3) can be written: 

^{E+\P\ cose), (B.4) 

where 6 is the angle between P and P;. Let us now take the limit of a massless tau, i.e., 
E = \P\, then we get: 

tE 

^(1 + cos^). (B.5) 

Assuming that the direction of the lepton is chosen isotropically^ in the rest frame 
of the tau, implies that the distribution of cos^ is flat. Prom that we can calculate the 
distribution of y = 1 + cos 9 as: 



fyiy) = o' (B-6) 



1 

2' 

where y £ [0, 2]. We arc interested in the fractional lepton energy with respect to that of 
the tau, in the Xi rest frame, i.e., the product z = where x by construction follows the 
distribution (A. 4) with i = 2,3. The distribution fziz) is now given by: 

„ , , 5-9^2+4^3 

fz{z) = . (B.7) 

With all decay products from the neutralino massless, we can simplify the invariant mass 
formula: 

M^b = (Pa + Pb)^ = 2EaEb{l - cos <^), (B.8) 

where (p is the angle between Pa and Pi,. If we then assume that particle 6 is a tau and 
decays to a charged lepton d, the resulting invariant mass Mad is given by: 

Mad = Mm (B.9) 
where q = -y/i and has a distribution fq{q) = ^{5 — 9q'^+4q^). Mad then has a distribution: 

fM^Maa) = C IMlMdMM^^^ (BIO) 



®Since the tau will in general be left-handed, this is not a very accurate assumption. For a discussion of 
the effects of tau polarization, see below. 
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Figure 21: Effects of cuts on invariant mass distributions involving taus. Left: No cuts. Right: 
Cuts as indicated by Eq. (B.ll). Long-dashed (green) curve: two tau jets are combined; sohd 
(blue): one directly produced lepton combined with a tau jet; short-dashed (blue): a tau together 
with a lepton from a tau decay. 

We see from Eq. (B.7) that there are many low-energy leptons coming from the tau decay. 
These low-energy leptons will, however, not pass the pT cuts we impose and therefore we 
need to take the cuts into account when calculating the distribution (this is not required 
for the direct decay to leptons since the amount of low energy leptons there is small) . The 
simplest thing to do to include some cuts, is to constrain the lower limit in the integral of 



The introduction of this cut does not preserve the normalization of the distribution; 
up to this point all distributions have been normalized to unity as one would expect for 
probability distributions. However, this is no longer true for Eq. (B.ll), moreover, the new 
normalization will depend on the distribution to which this kind of cut is applied, due to 
the larger impact on the lower end of the distribution. 

Since we are more concerned about the shape rather than the exact values from the 
distributions, there is no need to normalize Eq. (B.ll). Note also that introducing the cuts 
this way is simple and gives a reasonably good result, but it is very ad hoc, the specific 
value 0.4 is adopted for no other reason than the fact that it gives a good approximation 
to the Monte Carlo results. 

We can also extend the analysis to the case when both leptons come from decaying 
taus by multiplying Mad with q = ^/z once again. The result of this is shown as the long- 
dashed curves of Fig. 20 and the distributions of Eq. (B.IO) and Eq. (B.ll) are shown as 
the short-dashed curves. 

So far we have ignored any effects of tau polarization. However, since the neutralino 
will in general decay to left-handed taus, the spin of the tau will affect the mass distribution. 
More precisely, we get an extra factor ^(1 + cos 6) which amounts to replacing Eq. (B.6) 
with fy{y) = |. The result of this can be seen in the thin green curves of Fig. 20. The 



Eq. (B.IO): 




(B.ll) 
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Figure 22: Invariant mass distributions expected from the higher RPC chain. Left: No cuts. 
Right: Cuts as indicated by Eq. (B.ll). Sohd: the familiar triangular di-lepton distribution; short- 
dashed: the di-tau-jet invariant mass; long-dashed: lepton-tau-jet mass distribution, expected when 
one of the taus decays leptonically. 

tau polarization pushes the distributions to higher energy but the effect is smaller than 
the uncertainties regarding the effects of the cuts. Given that PYTHIA also does not 
take polarization into account, we use the distributions without polarization effects for the 
comparison with Monte Carlo. 

In the right panel of Fig. 20 we see the result of simulating the cuts with a constraint 
on the integral as done in Eq. (B.ll). As mentioned before, this suppresses the lower 
end of the distributions and changes the normalizations. We can see that the change in 
normalization is different without (blue curves) and with tau polarization effects (green 
curves). Including tau polarization effects, the distributions are pushed towards higher 
masses and are therefore less affected by the cuts. 

If the tau instead decays hadronically, we need to estimate the energy of the resulting 
tau-jet (rjet). The simplest thing to do is to assume a three-body decay of the tau where 
the neutrino gets a fraction c of the tau mass and the rest of the energy goes into the rjet- 
The distribution of c will then follow (A. 4) for i = 2,3. We get the invariant mass of a rjet 
combined with a lepton, Mt-.^j as: 



where Mab and c are as defined above. In line with the above procedure we introduce the 
cuts by requiring \/l — c > 0.4. 

This procedure can be repeated to produce the distribution of two rjet's as well as 
lepton- Tjet distributions where the lepton comes from the decay of a tau. The result of this 
is shown in Fig. 21. 

For comparison, in Fig. 22 we show the expected signals from the RPC chain. One 
clearly sees the similarities between the RPV and RPC signals in channels including rjet's. 
This makes the study of these channels much more difficult. 



Mr.^J = MabVl - C, 
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